Phonons in the /3-tin, Imma, and sh phases of Silicon from ab initio calculations 
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We present a new interpretation of measured Raman frequencies of a high-pressure structure of 
Silicon which was assigned previously to the /3-tin phase. Our results show that the /3-tin— >Imma— »sh 
phase transitions have been already indicated in this experiment which was performed before the 
discovery of the Imma phase. We have calculated phonon-dispersion curves for the /3-tin, Imma, 
and sh phases of silicon using the plane-wave pseudopotential approach to the density-functional 
theory and the density-functional perturbation theory within the local density approximation. With 
the new assignment, the calculated phonon frequencies display an excellent agreement with the 
experimental data, and can be also used to determine precisely the transition pressure for the 
Imma^(3-tin phase transition. The sh—>Imma transition is accompanied by soft modes. 

PACS numbers: 61.50.Ks 63.20.Dj 71.15.Mb 71.15.Nc 



In 1993 a new high-pressure phase was found between 
the /3-tin (Sill, body centered tetragonal structure) and 
the sh (SiV, simple hexagonal structure) phase of silicon 
0. This phase was called Imma phase after its space 
group and has a body-centered orthorhombic structure 
(SiXI). In the diffraction experiment, the phase transi- 
tions to and from the Imma phase appear as of first or- 
der with a discontinuity in the volume and the lattice pa- 
rameters • Previous ab initio calculations of the phase 
transitions /3-tin— >Imma— »sh have indicated both phase 
transitions to be of second order Q or both of first order 
0. A theoretical investigation of elastic stabilities has 
shown a second-order phase transition /3-tin— > Imma 5]. 
From group-theoretical arguments 0, the phase transi- 
tion /3-tin— > imma can be of second order, whereas the 
phase transition imma— >sh has to be of first order. This 
conclusion was strongly supported by an ab initio calcu- 
lation in our previous work [jj. 

An intermediate phase similar to the Imma phase, with 
a structure more general than the /3-tin and the sh phase 
was considered theoretically before || with a prediction 
that the phase transitions might be accompanied by soft 
phonon modes. Since the Imma structure results from 
a distortion and a relative sublattice shift of /3-tin along 
the c-direction, which corresponds to a soft T phonon 
displacement [9j , it is worthwhile to investigate the pres- 
sure dependence of the phonon-dispersion curves. Es- 
pecially, it was speculated that the superconductivity of 
these metallic high-pressure phases of silicon might be 
enhanced by these soft modes 0, 0] . Available experi- 
mental and theoretical data for superconducting proper- 
ties [II HQ EES 

did not consider the Imma phase, 
because its existence was not known at that time. In 
the experiment, the superconduc ting temperature as a 
function of pressure shows a kink [ll|, [ljj [l6| within the 
assumed stability range of the sh phase. In the corre- 
sponding theoretical work this kink was traced back to 
some soft modes from calculations of a few phonon modes 
under pressure for the sh structure with a frozen-phonon 



technique 13]. However, in our opinion this kink corre- 
sponds to the phase transition Imma— >sh. A first step 
towards the determination of the superconducting tran- 
sition temperature by calculating the electron-phonon 
coupling is an analysis of the phonon frequencies of all 
phases. 

In this contribution we examine the phonons of the 
high-pressure phases /3-tin, Imma, and sh and their be- 
haviour near the phase transitions. In particular, we 
are interested in soft modes and the pressure-dependent 
phonon frequencies. An indication of a soft-mode be- 
haviour has been already found experimentally: The low- 
energy optical phonon frequency of the /3-tin phase at 
the T-point decreases with increasing external pressure 
[l7j ]. Restricting ourselves to just the /3-tin phase we did 
not find a vanishing frequency even beyond the stabil- 
ity range of this phase |lg, LUj , but we find a lowering 
of the lower optical phonon frequency, whereas the de- 
generate upper ones increase. In the following we want 
to include especially the Imma and the sh phases in our 
consideration. 

This paper is organised as follows: First, we give a 
short overview of the methods on which our calculations 
are based. Second, we describe in brief the structures, 
the relaxation and the pressure dependence of the struc- 
tural parameters of the three phases. Third, phonon- 
dispersion curves of the three phases are presented as 
well as an analysis of the phonon frequencies at the zone- 
center and at a zone-boundary point as a function of the 
external pressure. Finally, we draw a conclusion. 

The calculations have been carried out using a plane- 
wave pseudopotential scheme (ijj within the density- 
functional theory EHl and the density- functional per- 
turbation theory |23J, |2J, |2jj( implemented in the PWscf 
package j2(| . The ion core of silicon has been described 
by a norm-conserving pseudopotential 1271 l28l created 
following a scheme by v. Barth and Car [29j, which is 
described in [^O]. The exchange-correlation energy is 
calulated with the use of the local-density approximation 
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FIG. 1: Relaxed lattice parameters of the BCO cell as a func- 
tion of the external pressure. The filled symbols mark the 
stability range of the Imma phase between the /3-tin phase 
(A = 0.25) and the sh phase (A = 0.5). The dashed line 
marks the ideal value b/c = \/3 for the SH structure. 



[3lL |32( . The pressure corresponding to a given structure 
is obtained by evaluating the stress tensor (3^, 01 ■ We 
have used a kinetic-energy cutoff E r „ t of 40 Ry and 1165 
special Monkhorst-Pack points [35j in the irreducible 
part of the Brillouin zone with a Methfessel-Paxton |3jj 
smearing width of 0.03 Ry, because of the metallic char- 
acter of the system. With these parameters the conver- 
gence is better than 0.01 mRy for energ y di fferences, and 
the pressure related to the Pulay stress [33,|38| is smaller 
than 1 kbar. The phonon frequencies have been calcu- 
lated on a discrete mesh of 18 points in the irreduzible 
Brillouin zone, and the phonon-dispersion curves have 
been calculated by Fourier interpolation. 

All these three structures have been determined using 
the same body-centered orthorhombic (bco) cell with lat- 
tice constants a ^ b ^ c and two atoms at (0, 0, 0) and 
(0,0. 5&, Ac) to which we refer as BCO. In this cell the 
/3-tin phase can be reproduced by choosing a = b and 
A = 0.25 and the structure of the sh phase by b = V3c 
and A = 0.5. The relaxation of the lattice parameters 
has been carried out by evaluating the total energy as a 
function of the lattice constants. The equlibrium lattice 
parameters for a fixed volume have been determined by 
the condition p x — Py — Pz where p x , p y , and p z are 
the negative diagonal components of the stress tensor. 
Therefore, the external pressure can be taken from the 
stress tensor since it is diagonal with equal components 
for a relaxed structure. For the BCO structure, the in- 
ternal parameter A has been relaxed in addition to the 
lattice constants, while for the /3-tin and the sh phases 
A was fixed at the corresponding value. In the latter 
case the ideal relation a = b has been reached for the 
/3-tin phase but not exactly for the sh one. On the other 
hand, the deviations from 6 = \/3c are negligible and do 
not affect the calculation of the phonon frequencies. The 
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FIG. 2: Phonon-dispersion curves for the /3-tin phase 
(V = 188 a|,98 kbar, dashed), the Imma phase (V = 
182 ag, 140 kbar, solid), and the sh phase (V = 
176 a|,183 kbar, dotted) at pressures within the stability 
range of each phase. The corresponding directions are from 
left to right: A [001], A [010], E [100], [111], and [011]. 



results corresponding to these structures are denoted as 
BCT and SH. 

In fact, for a free relaxation of the BCO cell we obtain 
the lattice parameters of the /3-tin phase for large vol- 
umes (small pressures), the ones of the sh phase for small 
volumes (large pressures) and between these phases the 
Imma phase with 0.25 < A < 0.5 as apparent in Fig. ^ 
Here, with A being discontinuous the phase transition /3- 
thi—tlmma seems to be of first order, while the Imma— *sh 
transition seems to be of second order because of a miss- 
ing discontinuity of A. The enthalpy barrier between the 
Imma and the sh phase is only a few meV Q, which is 
less than the precision of our calculations. Therefore, the 
determination of the order of the phase transitions from 
the behaviour of the lattice parameters as a function of 
the pressure or total energy is unreliable. Nevertheless, 
we can estimate the transition pressures from Fig.^ and 
we find that the /3-tin— > Imma transition occurs between 
115 and 118 kbar and the Jmmo^sh one between 142 
and 143 kbar, which is in good agreement with the ex- 
perimental values if those data are evaluated in the 
same way (for further details and a review of available 
experimental and theoretical results see 0). In the fol- 
lowing we will inspect the phonon spectra as a possible 
source for clearer conclusions. 

Phonon-dispersion curves have been calculated for the 
three phases at various pressures. For each structure at a 
given volume within the stability range of the correspond- 
ing phase the phonon frequencies along selected high- 
symmetry wave-vector directions are shown in Fig. [2] 
Compared to the curves of the /3-tin phase some degen- 
eracies of freqencies phase at high-symmetry points are 
lifted in the Imma because of the symmetry of the this 
phase being lower than that of /3-tin, see, e.g., at the T- 
and the X-point. For the sh phase new degeneracies ap- 
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pear because of the higher symmetry of this phase with 
respect to the Imma phase. Surprisingly, while for the 
/3-tin and the Imma phase the polarisations of the modes 
are very similar (even a simultaneous exchange of the po- 
larisation of two modes along some directions), some ma- 
jor changes appear between the Imma and the sh phase. 
For example, around the middle of the XUr-direction 
the third mode (counted from high to low frequency at 
U) has a longitudinal acoustic polarisation whereas it has 
a transverse optic polarisation in the case of Imma and 
/3-tin. 

Like in the cd phase of Si the high-energy optic 
frequencies increase and the low-energy acoutic ones de- 
crease with increasing pressure. In the present case, an 
exceptionally strong decrease is found at the V- and the 
S-points. In order to find some critical phonon softening 
we have investigated the frequencies as a function of ex- 
ternal pressure at these points. For this reason we were 
not only taking into account the fully relaxed BCO struc- 
ture but also the BCT and SH structures. The results 
are presented in Fig. EI Here the three phases are visible 
again, meaning that at small pressures the frequencies of 
BCO are identical with the ones from the BCT structure 
corresponding to the stability of the /3-tin phase. At high 
pressures the same correspondence can be stated for sh. 
The shift of the frequencies between /3-tin and Imma is 



much smaller than between Imma and sh which points to 
a second-order phase transition for /3-tin— >Imma and a 
first-order phase transition Imma^sh. We assume that 
the change will become more abrupt between Imma and 
sh for more accurate ground-state calculations which de- 
scribe the energy barrier between these phases with more 
precision. Thus, we expect a discontinuity of the fre- 
quency at this pressure from a more accurate calculation. 

In principle, the softening of a phonon mode indicates 
an instability against a displacive phase transition. We 
find a softening with decreasing pressure in the SH struc- 
ture for the sh phase at the zone-enter T-point and at the 
zone-boundary S-point around the sh—> imma phase tran- 
sition. The instabilities at both points are related to an 
instability at the L-point of the Brillouin zone of a simple- 
hexagonal cell using the corresponding monatomic struc- 
ture. A softening at a point on the zone boundary cor- 
responds to a doubling of the unit cell. In fact, in the 
ideal SH symmetry the BCO unit cell contains two sh 
unit cells. We can estimate the transition pressure for 
the sh—>Imma transition from these soft modes and find 
the values 143.9 and 146.3 kbar from the softening at the 
T- and the S-point, respectively. These transition pres- 
sures are in good agreement with the ones obtained from 
the relaxed lattice parameters. 
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FIG. 3: Pressure dependence of the phonon frequencies at the T- and the S-point for the BCT-, BCO-, and SH structure. 
Imaginary frequencies are drawn along the negative frequency- axis. The dots represent calculated values. Lines are guides to 
the eye in the case of BCO, and in the case of BCT and SH are derived from the linear function co 2 (p). Reference values from 
other calculations (TH) stem from 0] (square), (39|l (diamonds), and pfll (triangles). Experimental data from \vj\ are denoted 
as crosses. The pictures of the Brillouin zones (bco upper panel, sh lower panel) are used as an illustration. 



We find an excellent agreement also with the avail- 
able experimental data of Olijnyk 0] who measured 
the pressure-dependent phonon frequencies for the /3- 
tin structure at a time when the Imma phase was yet 
unknown. Wrongly assuming the /3-tin phase, theoret- 
ical work has not been able to reproduce these data 



[lil pil liof . Our results agree completely with those of 
Ref. [IT! in the low-pressure (/3-tin) and the high-pressure 
(sh) range and describe the slope in between (Imma). 
Therefore, Olijnyk was actually observing the phonon 
frequencies of all these three phases. 

In the Imma phase the TO modes are split into two 
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branches which are both Raman active. Since the peak 
in the experimental Raman spectrum is very broad and 
the frequency of the maximum is just a little bit lower 
than our calculated value we speculate that the second 
peak has a lower intensity and is unresolved from the 
high-frequency one. Besides, there might be a second- 
oder background in the spectra: If the overtone spectra 
are dominant in the /3-tin, Imma, and sh phases as they 
are in the diamond phase, then the spectra correlate with 
the DOS of Fig. if the frequency axis is stretched by 
a factor of 2. In fact, the experimental intensity near 
500 cm -1 has a strong positive slope and is larger in the 
/3-tin than in the Imma or sh phases. 

In conclusion, we have investigated the phonon dis- 
persion curves for the /3-tin, Imma, and sh phases of 
Silicon. While the order of a phase transition often is 
derived from the behaviour of the structural parameters 
or of the volume, the order of the pressure-induced (3- 
tia-^Imma— >sh phase transitions can also be obtained 
from the pressure dependence of the phonon frequencies. 
Here we find a small shift of the frequencies between the 
/3-tin and the Imma phase and a large one between Imma 
and sh pointing to a second-order and first-order phase 
transition, respectively Our calculations are able to re- 
produce the experimental data of Olijnyk ^3] perfectly. 
This results in the insight that in the Raman spectra the 
Imma and the sh phases are visible rather than the as- 
sumed /3-tin phase. There is a splitting of the Raman 
mode at the /3-tin— >Imma transition as well as a critical 
softening of the LO mode at V and S at the sh—> Imma 
transition which has escaped experimental observation 
so far. Since the pressure-dependent behavior of these 
phases is known now, it should be a challenge to find 
the lower TO mode of the Imma and the sh phases and 
also the softening of the LO mode in the sh structure 
experimentally. 
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